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ABSTRACT 

Using an extensive sample of galaxies selected from the Sloan Digital Sky Survey Data Release 5, we com- 
pare the angular momentum distribution of active galactic nucleus (AGN) with non-AGN hosting late-type 
galaxies. To this end we characterize galactic spin through the dimensionless angular momentum parameter 
A, which we estimate approximately through simple dynamical considerations. Using a volume limited sam- 
ple, we find a considerable difference when comparing the empirical distributions of A for AGNs and non-AGN 
galaxies, the AGNs showing typically low A values and associated dispersions, while non-AGNs present higher 
A values and a broader distribution. A more striking difference is found when looking at A distributions in thin 
Mr cuts, while the spin of non-AGN galaxies presents an anti-correlation with Mr, with bright (massive) galax- 
ies having low spins, AGN host galaxies present uniform values of A at all magnitudes, a behavior probably 
imposed by the fact that most late-type AGN galaxies present a narrow range in color, with a typical constant A 
value. We also find that the fraction of AGN hosting galaxies in our sample strongly depends on galactic spin, 
increasing dramatically for decreasing A. For AGN host galaxies we compute the mass of their supermassive 
black holes and find that this value tends to be higher for low spin galaxies, even at fixed luminosity, a result 
that could account, to a certain extent, for the spread on the luminosity black-hole mass relation. 
Subject headings: galaxies: active — galaxies: fundamental parameters — galaxies: spiral — galaxies: statis- 
tics — galaxies: structure 



1. INTRODUCTION 

Studies of central gas and stellar kinematics of nearby 
galaxies have revealed the presence of supermassive black 
holes (SMBHs) in almost all massive galaxies. According 
to the most accepted paradigm (Rees 1984), galactic nuclear 
activity is directly related with the fueling process of these 
SMBHs, which strongly correlates to the overall structure 
and evolution of the galaxies. For instance, active galactic 
nucleus (AGN) activity seems to be triggered only in galax- 
ies with massive bulges and intermediate colors, as pointed 
out by Choi et al. (2009), henceforth CWP09. These au- 
thors show how the necessary conditions for AGN activity 
depend on galaxy morphology; for early-type galaxies low ve- 
locity dispersion and blue colors are preferred, and indicative 
of available gas supply for the central black hole, while in- 
termediate velocity dispersion, intermediate-color, and more 
concentrated late-type galaxies are more likely to host AGNs, 
reinforcing the idea that AGN host galaxies are intermediate 
objects between the red sequence and the blue cloud (Kauff- 
mann et al. 2003). 

If most galaxies present SMBHs in their centers, a funda- 
mental question is why not all of them host AGNs when gas 
supply is available. This is a question with no simple answer 
given the complexity of the processes involved in the forma- 
tion and evolution of galaxies, but important clues have been 
proposed to elucidate the mechanisms that trigger AGN activ- 
ity (Martini 2004). Major mergers of gas-rich progenitors and 
galaxy interactions are two commonly invoked mechanisms 
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that can provide gas to the center of the galaxies and fuel dor- 
mant SMBHs (Hopkins et al. 2006; Ellison et al. 2008). Tidal 
torques exerted by non-axisymmetric features such as large- 
scale and nuclear bars (Jogee 2006), recycled gas from dying 
stars in the inner kiloparsecs of the galaxy (Ciotti & Ostriker 
2007), and enhanced viscosity of the turbulence associated 
with supernova explosions (Chen et al. 2009) are just some 
other mechanisms that can effectively trigger nuclear activity. 

A common requirement for any of the mentioned mech- 
anisms is the removal of specific angular momentum in or- 
der for the material to collapse to the center, and operate as 
fuel. Therefore, a relation between the overall galactic A spin 
parameter and the onset of AGN activity is expected, with a 
higher frequency of AGNs in low spinning galaxies. Statisti- 
cally, AGN hosting galaxies should form a low A population, 
in comparison to normal galaxies. Low angular momentum 
is expected not only in the growing phase of SMBHs, but has 
also been invoked to explain the formation of massive black 
hole seeds themselves. In this scenario (Eisenstein & Loeb 
1995; Koushiappas et al. 2004), the direct collapse of low 
angular momentum material results in a typical inverse pro- 
portionality between the SMBH mass and the total angular 
momentum (Haehnelt et al. 1998; Colgate et al. 2003). 

In this work, we study the spin distribution of AGN host 
galaxies in contrast with non-AGNs, parameterizing the an- 
gular momentum through the dimensionless spin parameter 



where E, M and L are the total energy, mass and angular 
momentum of the configuration, respectively (Peebles 1969). 
Our aim is to determine empirically if the AGN population 
effectively presents a low spin distribution when compared 
with non-AGNs and corroborate the results arising from the- 
oretical studies concerning the relation between A and the 
SMBH mass. To this end, we use a sample of galaxies se- 
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lected from the Sloan Digital Sky Survey (SDSS) and apply 
a simple model to obtain an empirical estimation of A. The 
large volume limited sample allows a thorough coverage of 
parameter space, which permits an empirical determination 
of the leading galactic physical parameters driving the AGN 
phenomena and the observed scalings with SMBH mass. The 
physical model is presented in § 2, and the sample selection 
criteria in § 3. In § 4 we present our general results and the 
final conclusions are summarized in § 5. 

2. ESTIMATION OF THE SPIN FROM OBSERVABLE 
PARAMETERS 

In Hernandez & Cervantes-Sodi (2006), we derived a sim- 
ple estimate of A for disk galaxies in terms of observational 
parameters and showed some clear correlations between this 
parameter and type defining structural parameters. Here we 
briefly recall the main ingredients of the simple model. The 
model considers only two components, a disk for the baryonic 
component with an exponential surface mass density S(r), 



(2) 



where r is a radial coordinate, and Eq and are two con- 
stants which are allowed to vary from galaxy to galaxy, and 
a dark matter halo having a singular isothermal density pro- 
file p{r), responsible for establishing a rigorously flat rotation 
curve Vd throughout the entire galaxy: 



p(r): 
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In this model we are further assuming that (1) the specific 
angular momentum of the disk and halo are equal (e.g. Mo et 
al. 1998; Zavala et al. 2008); (2) the total energy is dominated 
by that of the halo which is a virialized gravitational structure; 
and (3) the disk mass is a constant fraction of the halo mass, 
F = Mi/ Mil. These assumptions allow us to express A as 



A: 



GMh 



(4) 



Finally, we introduce a baryonic Tully-Fisher (TF) relation: 
Md = ATpVd ^ (e.g., Gurovich et al. 2004). Taking the Milky 
Way as a representative example, we evaluate F and Ajf to 
obtain 



A = 21.8 



^d/kpc 



(yd/kms-i)^^/2- 



(5) 



We proved the accuracy in our estimation of A in Cervantes- 
Sodi et al. (2008), comparing the estimation using Equation 
(|5]l to values arising from numerical simulations of six dis- 
tinct groups, where the actual value of A is known, as it is 
one of the parameters of the simulated galaxies, and where 
this can also be estimated through Equation (|5]l, as baryonic 
disk scale lengths and disk rotation velocities are part of the 
output. The result was a one-to-one correlation, with very 
small dispersion and no bias, leading to typical errors < 30%, 
which we include throughout. The use of this estimate for 
large SDSS samples led us in Hernandez et al. (2007) and 
Cervantes-Sodi et al. (2008) to a first empirical derivation of 
the distribution of A parameters in the local universe and a 
measure of the scalings of this parameter with local environ- 
ment and halo mass. These results were latter confirmed by 
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Fig. 1 . — Galaxy distribution in tlie BPT diagram, pure AGN galaxies lying 
above the solid line (Kewley et al. 2006), star-forming galaxies below the 
dashed line (Kauffmann et al. 2003), and composite objects in between. 

Berta et al. (2008) using an estimate of A through Equation 
(|5]l, also used by Gogarten et al. (2010) as an unbiased first- 
order A measurement. 

3. THE SDSS SAMPLE 

The sample of galaxies used in this work is an enhanced 
version of the study by CWP09. The galaxies come from the 
large-scale structure sample, DR4plus, using the New York 
University Value-Added Galaxy Catalogue (VAGC; Blanton 
et al. 2005), which is a galaxy subset of the SDSS Data 
Release 5 (Adelman-McCarthy et al. 2006), and the Korea 
Institute for Advanced Study-VAGC (Choi et al. 2010). A 
detailed description of the sample can be found in Choi et 
al. (2007) and CWP09. The rest-frame absolute magnitudes 
of individual galaxies were computed in the r band, using 
Galaxy reddening correction (Schlegel et al. 1998) and K- 
corrections as described by Blanton et al. (2003), with a mean 
evolution correction given by Tegmark et al. (2004). We 
define a volume-limited sample with a lower magnitude cut- 
off M,. = -19.0, with redshifts 0.020 < z < 0.086, and a total 
number of 87,590 galaxies. Disk scale lengths are estimated 
from one-component light fits, which implicitly account for 
the presence of low angular momentum material in bulges. 

In order to compare the spin of AGN and non-AGN galax- 
ies, we use only Type II AGNs to avoid complications due to 
the presence of a bright nucleus. The AGN galaxies were seg- 
regated from star forming galaxies using the flux ratio BPT 
diagram of Balmer and ionization lines (Baldwin et al. 1981). 
The AGNs were selected among galaxies fulfilling the defini- 
tion by Kewley et al. (2006): 

0.61/(log([Nn]/Ha)-0.47)+1.19<log([Oni]/H/3), 
0.72/(log([Su]/Ha)-0.32)+1.30<log([Om]/H/3), 
0.73/(log([Oi]/Ha) + 0.59)+1.33<log([Om]/H/3), 

while star-forming galaxies were selected using the selec- 
tion criteria by Kauffmann et al. (2003): 

0.61/(log([N„]/Ha)-0.05)+1.30<log([Om]/H/3), 
0.72/(log([S„]/Ha)-0.32)+1.30<log([Om]/H/3), 
0.73/(log([Oi] /Ha) + 0.59) + 1 .33 < log([Om] /H/3), 
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Fig. 2. — (a) P(A) distribution for late-type galaxies in the sample segregated in normal, composite and AGN galaxies, (b) P(A) distribution for normal galaxies 
split according to absolute magnitude, (c) P(A) distribution for AGN galaxies split according to absolute magnitude, (d) Median A value of AGN galaxies and 
the control sample as a function of Mi. 



requiring a signal-to-noise ratio > 6. The galaxies resid- 
ing in the area between the previous demarcation lines (solid 
and dotted lines, respectively, in Figure [1} are identified as 
composite objects, which contain AGN as well as extended 
Hn regions. Our non-AGN sample consist of those galaxies 
excluding AGNs and composite galaxies, without signal-to- 
noise cut. 

To estimate the black hole mass (Mbh), we use the Mbh- 
stellar velocity dispersion relation by Tremaine et al. (2002): 

log(MBH/Mo) = 8.13 + 4.02 x log(cr/200 km s"'), (6) 

where a simple aperture correction to the stellar velocity 
dispersion was applied. Since velocity dispersion measure- 
ments lower than the instrumental resolution are not robust, 
we will focus on AGNs with Mbh > lO^^M©. 

In order to discriminate between elliptical and disk galax- 
ies, we used the prescription of Park & Choi (2005) in which 
early (ellipticals and lenticulars) and late (spirals) types are 
segregated in a m - r color versus g- i color gradient space and 
in the concentration index space. They extensively tested the 
selection criteria through a direct comparison of visually as- 
signed types for a large sample of several thousand galaxies. 
The specific selection criteria can be found in Park & Choi 
(2005). For more details the reader is referred to CPW09. 

4. RESULTS 



Our Equation (|5]l requires the rotational velocity to calcu- 
late A, which is inferred from the absolute magnitude in the 
r-band and using a TF relation (Pizagno et al. 2007). To 
avoid the problem of internal absorption in edge-on galaxies 
(Cho & Park 2009), and consequently underestimating rota- 
tional velocities, we limit the sample to spiral galaxies having 
seeing-corrected isophotal axis ratios b/a> 0.6. 

Once the galaxies in our sample are segregated into early- 
and late-type galaxies, and into different activity classes such 
as normal, composite, and AGN hosting galaxies, we com- 
puted A using Equation (|5]) for all the late-type galaxies and 
obtained the spin distribution for the three different sub- 
classes: 41,662 normal, 2,273 composite and 1,026 AGNs. 
Figure 2(a) shows the distribution of A values for normal, 
composite, and AGN late-type host galaxies. We can see a 
clear difference between the normal and the AGN population, 
the former showing typically higher spin values with a large 
dispersion compared with the latter, which is a more coherent 
population of low spin galaxies. The distribution of galaxies 
classified as composite follows almost exactly the same dis- 
tribution as AGN galaxies, probably due to the low number of 
misclassified star-forming galaxies. Given that these galaxies 
are actually a mixed group of AGN and non-AGN galaxies 
we will not consider them for the rest of the analysis, and we 
will focus only in the difference between AGNs and normal 
galaxies. 

Theoretical (Shaw et al. 2006), as well as empirical (Her- 
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nandez et al. 2007) distributions of A, are commonly de- 
scribed by a log-normal function, in our case, the values de- 
scribing the non-AGN and AGN distributions are, respec- 
tively, Ao = 0.059 ± 0.006, ax = 0.400 ± 0.005 and Ao = 
0.042 ± 0.004, ax = 0.350 ± 0.003, showing a clear difference 
between the two populations. 

We use the two sample Komogorov-Smimov test to com- 
pare the distributions of normal and AGN host galaxies to 
determine if they are drawn from the same parent distribu- 
tion. The result of the test yields a negligible probability of 
P ^ X 10""**' that the two samples are drawn form the same 
parent distribution. 

If we look at the A distribution of normal galaxies, splitting 
the sample according to the absolute magnitude (Figure 2(b)), 
we clearly see that the distributions get tighter and present 
lower mean A values for brighter galaxies, the same behavior 
we reported in Cervantes-Sodi et al. (2008), where we found 
that massive galaxies tend to have low spin dark matter halos 
and show low dispersion, while low mass galaxies have higher 
A values and high dispersions around the median value, a re- 
sult later confirmed by Berta et al. (2008). A striking differ- 
ence appears for the case of AGN host galaxies (Figure 2(c)), 
showing low mean A values and tight distributions for all Mr- 
This result indicates that only low spinning galaxies, at all 
mass ranges, are suitable AGN hosts in the local universe. 

Previous work has already pointed out that the recurrence of 
AGN activity is a function of internal properties such as abso- 
lute magnitude, color, and velocity dispersion. Kauffmann et 
al. (2003) state that massive galaxies with old stellar popula- 
tions frequently host AGNs and CWP09 found that the frac- 
tion of AGN increases with galaxy luminosity, and particu- 
larly for late-type galaxies, this fraction increases for interme- 
diate color systems, with M-r = 2~2.4 and intermediate ve- 
locity dispersion (~ 130 km s"'). In this context, most of our 
AGN hosting galaxies are restricted to having intermediate 
colors, a restriction that also constrains the value of A given 
that color anticorrelates with A (Hernandez & Cervantes-Sodi 
2006). To make a fair comparison between AGN and non- 
AGN galaxies, we constructed a control sample by selecting 
non-AGN galaxies which match one to one the morphology, 
absolute magnitude, color and concentration index of each 
galaxy in the AGN sample, guaranteeing that the distributions 
of these parameters are the same for both populations. Figure 
2(d) shows how the typical A value of AGN hosting galaxies 
tends to be lower than the one presented by the control sam- 
ple at all M, values. The uncertainties (throughout this paper) 
represent 1 a confidence intervals determined by the bootstrap 
resampling method. 

A natural question arising from our previous results is 
whether or not the AGN fraction depends on galactic spin. 
Previous studies have shown the dependence of the AGN frac- 
tion on the brightness of the parent galaxies (Mercurio et 
al. 2010, CWP09), or their associated stellar or total mass 
(Pasquali et al. 2009). Figure 3(a) shows the fraction of 
AGN-Hcomposite (/AGN+Composite) galaxies as a function of ab- 
solute magnitude in different A cuts, where a moderate in- 
crease of /AGN+Composite IS present for increasing brightness, 
while a large difference in this fraction is evident for the dis- 
tinct A cuts. Figure 3(b) gives a complementary study for 
/AGN+Composite as a function of galactic spin for thin M,. slices, 
where the dependence on A appears stronger than the case of 
M,-, with all curves consistent with a single underlying depen- 
dence to within errors, for all magnitude bins. This shows 
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Fig. 3. — AGN fraction as a function of(a) absolute magnitude with cuts of 
A and as a function of A with cuts in magnitude (b) and color (c). 



that the driving physical ingredient is the value of A rather 
than magnitude. 

CPW09 have already shown that late-type galaxies present 
weak dependence of /ACN+Composite on magnitude or mass (ac- 
counted for through velocity dispersion) in addition to m - r 
color, and in the preceding section we found that the value 
of A is intimately linked to the color of the galaxy in ques- 
tion. The /AGN+Composite dependency on A for narrow u-r bins 
is presented in Figure 3(c), where even at fixed color, a sub- 
stantial decrease of the AGN fraction results as A increases, 
at least in the range of colors where the highest occurrence of 
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AGNs is found. We see that regardless of mass or color, high 
A systems will only very rarely present an AGN. 

The tight correlations of Mbh with host galaxy properties, 
such as the well-known one with the bulge luminosity (Ko- 
rmendy & Richstone 1995) or the tighter Mbh versus bulge 
velocity dispersion relation (Graham et al. 2011), have now 
been extended to more global properties, not related to the 
bulge component alone. Correlations with the spiral arm pitch 
angle of disk galaxies (Seigar et al. 2008), the mass of the 
dark matter halo (Ferrerase 2002; Booth & Schaye 2010), and 
the total gravitational mass of the host galaxies (Bandara et 
al. 2009) show how the correlations are not restricted only to 
the bulge component and point to a SMBH-bulge-total-mass 
interrelation (Volonteri et al. 201 1). Given the important role 
the spin plays in establishing the morphology and present-day 
structure of disk galaxies, and our previous result that shows 
the direct influence of A on determining nuclear activity, it 
is logical to expect a correlation between A and Mbh for the 
active galaxies in our sample; more specifically, an anticorre- 
lation. 

Let us consider the simple criterion by Larson (2010) for 
the formation of black holes, where gas in a forming galaxy 
with a column density above a critical value Scrit, goes into 
forming a central black hole. If the gas in the forming galaxy 
is distributed radially like the total mass with an isothermal 
density profile given by Equation (|3]l, with a gas to total mass 
ratio /g, the surface density of gas at any radius r in the disk 
will be Sg(r) = f^V]^ /lirGr. For a rotationally supported disk, 
the scale radius is given by ~ Arvir, where ryir is the virial 
radius of the system. This allow us to express the mass of gas 
within the critical column density criterion as 



My(S < I],rit) = 



/,m2 



7rSci-it(Arvir) 



(7) 



This expression, that does not attempt to be more than a di- 
mensional analysis estimate, implies that at fixed total mass 
the mass of gas that could end up forming a black hole is 
Mbh oc A"", with n = 2. Theoretical studies report different 
values for n that fall in the range 1-4 (Colgate et al. 2003; 
Koushiappas et al. 2004), all of them establishing an an- 
ticorrelation between the mass of the black hole and halo 
spin. Figure 4 (top panel) shows the empirical dependence 
of the inferred Mbh with the spin parameter for our sample 
of AGN host galaxies, showing the expected anticorrelation 
with A, in this case with a best fit described by log(MBH) = 
(-0.716 ± 0. 1 14) log(A) +(6.455 ± 0.155). In Figure 4 (bot- 
tom panel) we present the result of splitting the sample ac- 
cording to the absolute magnitude of the galaxies, where we 
can appreciate how at fixed A brighter (more massive) galax- 
ies host more massive black holes, with the anticorrelation 
with A present at all fixed Mr cuts. In this case, the leading 
effect is the total mass, with A introducing only a small but 
well-defined trend in the expected sense. 

5. CONCLUSIONS 

Using an extensive sample of galaxies extracted from the 
SDSS DR5, we show that the empirical distributions of the A 
spin parameter of AGNs and non-AGN host galaxies are qual- 



itatively and quantitatively different, the first distribution hav- 
ing lower A values and smaller dispersion around the median 
when compared with the latter Another striking difference 
between the two populations is that while the spin of normal 
galaxies presents a dependence on absolute magnitude, with 
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Fig. 4. — Upper panel: Mbh versus A for AGNs hosted by late-type galax- 
ies, showing the median value for each bin. Lower panel: the same as in the 
upper panel but showing the median Mbh values splitting the sample accord- 
ing to the absolute magnitude. 

a systematic decrease of A with increasing luminosity, AGN 
host galaxies present at all M,- low A values. This result, in 
addition to the increase of the AGN fraction for decreasing 
A, highlights the requirement for galaxies to have low spin in 
order to host AGNs and should be taken into account in semi- 
analytic models that exclusively use the halo or stellar mass 
to reproduce the fraction of galaxies belonging to different ac- 
tivity classes (e.g., Wang & Kauffmann 2008; Fontanot et al. 
201 1). Seeding AGNs into modeled dark matter halos should 
primarily consider A, in order to better model the observed 
universe. 

For the AGN sample we found that the inferred mass of the 
SMBHs shows a weak dependence on the spin of the host- 
ing galaxies, with increasing mass for decreasing spin, a log- 
ical result if we need the gas to be accreted onto the central 
massive object, a process that can be more effectively accom- 
plished if the raw material originally has low angular momen- 
tum. We recovered the leading Mbh-A relation at fixed abso- 
lute magnitude, with typically higher SMBH mass for brighter 
galaxies at fixed A, showing a double dependence of Mbh on 
Mr and A. Given that the trend with the spin is present at all 
M, bins, this could account, to a certain degree, for the disper- 
sion in the well established MBH-luminosity relation. 
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by grant from the Kyung Hee University (KHU-20100179) 
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